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Abstract 

Aeroelastic stability analyses were performed to insure 
structural integrity of two counterrotating propfan blade 
designs for a NAVY / Air Force / NASA cruise missile 
model wind tunnel test. This analysis predicted if the 
propfan designs would be flutter free at the operating 
conditions of the wind tunnel test. Calculated stability 
results are presented for the two blade designs with 
rotational speed and freestream Mach number as 
parameters. An aeroelastic analysis code ASTROP2 
(Aeroelastic STabilitv and Response Qf Propulsion 
Systems - 2 Dimensional Analysis), developed at LeRC, 
was used in this project. The aeroelastic analysis is a modal 
method and uses the combination of a finite element 
structural model and two dimensional steady and 
unsteady cascade aerodynamic models. This code was 
developed to analyze single rotation propfans but was 
modified and applied to counterrotation propfans for the 
present work. Modifications were made to transform the 
geometry and rotation of the aft rotor to the same 
reference frame as the forward rotor, to input a non- 
uniform inflow into the rotor being analyzed, and to 
automatically converge to the least stable aeroelastic mode. 

Nomenclature 

[A] generalized aerodynamic matrix 

A 1 ,..,A m torsion contributions of the m normal modes 
about the reference axis 
(F(t)) nonaerodynamic nodal force vector 

{F q } time-independent nonaerodynamic force vector 

(G(t)] generalized motion-independent force vector 
h bending deflection 

i VT 

j mode index, j = 1, 2, 3,..., m 

[K^] generalized stiffness matrix 

[K 1 centrifugal softening matrix in physical 

coordinates 

[K({u))] nonlinear stiffness matrix in physical 

coordinates 

aerodynamic coefficients defined in Ref. 1 
L aerodynamic lift 

LCOEF aerodynamic coefficient matrix 
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m 

number of normal modes 

M a 

aerodynamic moment 

[M\ 

physical mass matrix 

[M ] 

g 

generalized mass matrix 

N 

number of blades 

[P] 

stiffness matrix defined in Eq. (14) 

{P({u),t)} 

aerodynamic nodal force vector 

{P({u 0 })} 

steady-state aerodynamic nodal force vector 

(q) 

vector of generalized coordinates 


amplitude of motion in generalized coordinates 

r 

interblade phase angle index, r = 1, 2, , N 

s 

blade length along the reference line 

t 

time 

(u) 

vector of blade deflections at grid points 

{u o } 

vector of steady state deflections at grid points 

[W] 

modal matrix expressed in terms of individual 
bending and torsion contributions along the - 
reference line, see Eq. (A3) 

w lv .,w m 

bending contributions of the m normal modes 
about the reference axis 

a 

torsional deflection 

y 

eigenvalue defined in Eq. (14) 

(AF(t)) 

perturbation nonaerodynamic nodal force 
vector 

(AP w (t)} 

perturbation aerodynamic nodal force vector 
(motion-independent) 

{AP«Au},t)} 

perturbation aerodynamic nodal force vector 
(motion-dependent) 

{Au(t)J 

vector of vibratory deflections at grid points 
measured relative to {u q ] 


structural damping ratio in j-th mode 

H 

real part of i (proportional to damping) 

V 

imaginary part of i *fy (flutter frequency) 


interblade phase angle 

toi 

modal matrix 

© 

frequency 

“5 

frequency of j-th mode 

Subscripts 

0 

> 

steady state value 

g 

generalized (modal) 

F 

values at flutter point 
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Superscripts 

differentiation with respect to time 
T transpose 

Introduction 

Propfan design and technology have sparked interest 
in military and civil applications due to the potential for 
very high fuel efficiency at cruise speeds up to Mach 0.85. 
Various studies have been conducted by NASA to evaluate 
propfan propulsion systems (Refs. 1-3). One can expect 
from these studies that a counterrotating propfan 
propulser for cruise missiles offers the potential for 
significant increases in range or payload compared to 
current turbojet or turbofan designs. The long-range 
conventional standoff weapon (LRCSW) program was 
established as a joint U.S. Navy/ Air Force program to 
develop an advanced, long-range cruise missile powered 
by a propfan engine. As a part of this program, a joint 
Navy/Air Force/NASA Propfan Missile Interactions Project 
was initiated to determine the effects of a propfan engine 
on missile aerodynamics in a wind tunnel test (Ref. 4). 
NASA Lewis Research Center (LeRC) designed and 
fabricated two propfan blade sets for this wind tunnel test. 
The blade sets were chosen to bracket the various designs 
proposed by industry for the full-scale propfan missile 
engines. The major test objective was to investigate 
structural and aerodynamic interactions between the 
missile and the propfans. The LeRC team objectives of the 
program were to design and fabricate two composite 
counterrotating (CR) propfan blade models which would 
have high aerodynamic efficiency and structural integrity 
when mounted on a generic cruise missile model (Ref. 5). 
The models were to be tested in the NASA Ames Research 
Center 14 foot wind tunnel. A sketch of the wind tunnel 
model is shown in Figure 1. Each hub has six blades. 

The aerodynamic and acoustic requirements of 
propfans have resulted in designs with thin, swept, and 
twisted blades of low aspect ratio and high solidity 
compared to conventional propellers. Experimental and 
analytical research has shown that these characteristics of 
propfans give them the potential to flutter in their 
operating range (Refs. 1-3). Past theoretical and 
experimental propfan investigations at LeRC have resulted 
in various computer codes used to predict aeroelastic 
stability and response. One such research computer code is 
ASTROP2 (Aeroelastic ST abilitv and Response <2f 
propulsion Systems - 2 Dimensional Analysis). This code 
was developed to analyze the aeroelastic stability of single 
rotation propfan blades (Ref. 1). In the present work, it 
was modified and applied to analyze the counterrotation 
propfan blade models that were designed for this project. 


The objectives of this paper are: (1) to document the 
results of the aeroelastic stability (flutter) analysis of two 
CR propfan blade designs developed at LeRC for a cruise 
missile model, and (2) to describe modifications made to 
the ASTROP2 code that allowed analysis of the CR rotors 
and resulted in increased automation and accuracy. 

Analytical Formulation 

This formulation is described in detail in Refs. 1 and 6. 
A summary is presented here. The aeroelastic equation of 
motion of a blade can be written as 

[Ml {ii) + [KJ {uj + [K({u))l {u} = (P({ul,t)} + (F(t)} (1) 

where (u) represents the blade deflections at the grid 
points, [Ml is the mass matrix, [KJ the centrifugal softening 
matrix, [K([u])] the nonlinear stiffness matrix, {P([u],t)} the 
equivalent aerodynamic nodal force vector, and (F(t)) the 
equivalent nonaerodynamic force vector. Because of the 
large steady deflections and the consequent need for the 
geometric nonlinear theory of elasticity in which the strain 
and displacement relations are nonlinear, the stiffness 
matrix [K({u})] is a function of nodal displacements and, 
hence, is nonlinear. This provides the additional geometric 
differential stiffness due to centrifugal stiffening terms. 
The displacement dependent centrifugal softening terms 
are included in the matrix [KJ. The rotation also 
introduces Coriolis forces, but these have been shown to be 
negligible for thin propfan blades in Ref. 7. Hence, they are 
not included in the present formulation. 

Linearization of Equations 

Equation (1) is generally nonlinear and is valid for 
calculating performance, stalled and unstalled flutter, 
forced response, steady-state deflections, frequencies and 
mode shapes. One solution method is to directly integrate 
in the time domain, but it is computationally inefficient. 
Common practice is to perturb the displacements about a 
steady- state configuration by writing, 

(u(t)} = (uq) + (Au(t)J 

(P({u},t)) = (P({u 0 ))} + (AP( (Auj,t) + {AP w (t)J (2) 

{F(t)J = [F 0 } + [AF(t)J 

where {u 0 ), {P({u 0 })), and {F<>} are the steady state values of 
[u], {P((u},t)J, and [F(t)) } respectively. The quantities 
(Au(t)), (AP({Au},t)), and {AF(t)J are perturbations from f , 

{ P( { u 0 } )} , and {F 0 ), respectively. The perturbed 
aerodynamic force is split into motion-dependent 
(AP({Au},t)} and motion-independent (AP w (t)) parts for 
convenience. Substituting (u(t)}, (P((u),t)), and (F(t)] from 
Eq. (2) into the nonlinear Eq. (1) yields two sets of 
equations: one steady-state equation for (u 0 l and another 
for the perturbation variable (Au(t)]. These are 
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a r = 2kc / N 


r = 0, 1,2,..., N-l 


( 11 ) 


[[KJ + [K({uo))]] {u 0 } = {PC(uo})} + {F 0 1 (3) 

[M] [Aii(t)) + [[Kg] + [K([uo ))]] (Au(t)) 

= (AP([Au],t)} + (AP w (t)} + (AF(t)} (4) 

The steady -state configuration for a given rotational 
speed and Mach number is obtained by solving the 
nonlinear Eq. (3). The stiffness matrix [K({u 0 ])l includes 
elastic stiffness and differential stiffness due to centrifugal 
stiffening loads and steady-state aerodynamic loads. Once 
the steady-state deflection and the effective total stiffness 
are known from Eq. (3), the natural frequencies and mode 
shapes are calculated by solving 

[M] [Au(t)} + [[Kg] + [K({u 0 })]] {Au(t» = 0 (5) 

which leads to the generalized mass matrix [M g ], the modal 
matrix [0], and the modal frequencies C0j. 

Flutter and Forced Response by the Modal Method 

The general vibratory motion can be expressed as a 
superposition of the contributions of the various normal 
modes: 

(Au(t)} * [0] {q(t)J (6) 

Substituting Eq. (6) into Eq. (4) and post-multiplying the 
result by [bF leads to 

[M g ] ( q }+ [Kg] Jq) = [A] (q) + (G(t)} (7) 

where 

[M g ] = [0] T [M] [0] 

Kgjj = Mgjj coj 2 (1 + 2iCjj) 

[A] {q} = [0] T (AP({Au),t)) (8) 

(G{t)) = [<J)] T [{AF(t)} + {AP w (t)}] 

The order of Eq. (7) depends on the number of modes 
included in Eq. (6). This number is determined by 
performing numerical experiments. The structural 
damping in each mode is introduced through the damping 
ratio Q where the modal index is j = 1,2,3, m. 

Assuming that the motion is simple harmonic at flutter, 
the generalized coordinate vector of the blade can be 
written as 

[q] = (qj (9) 

For a tuned rotor, in which all the blades are identical, 
the aeroelastic modes consist of the N individual blades 
vibrating with equal amplitudes with a fixed interblade 
phase angle between adjacent blades (Ref, 8). For the s th 
blade vibrating in the r th interblade phase angle mode, Eq. 
(9) can be written as 

(q) = {qj e^^) (10) 

where the phase angle between adjacent blades is given as 


In Eq. (8), the generalized aerodynamic matrix is 
represented by [A] for each interblade phase angle for 
simple harmonic motion of the blade. This matrix is 
defined in Appendix A. The motion-independent 
aerodynamic and nonaerodynamic forces are represented 
by {GK t)} for each interblade phase angle. These forces are 
included for completeness in the formulation, but no results 
on forced response will be presented in this paper. After 
setting (G(t)} = [0] and substituting Eq. (10), Eq. (7) becomes 

- CD 2 [M g ] [q 0 )+ [Kg] (q 0 ) = [A] (q 0 ) (12) 

Simplifying further, the flutter eigenvalue problem can be 
written as 

[P] (q^ = y[Mg] [q] (13) 

where 

[P] = [Kg] - [A] (14) 

i ify = ico = p + iv 

Flutter occurs when p > 0. co is the flutter frequency and yis 
the eigenvalue obtained from the solution of Eq. (13). 

ASTROP2 Computer Code 

A detailed description of the ASTROP2 code can be 
found in Refs. 1 and 6. The ASTROP2 code is based on a 
fully coupled normal mode structural model and 2D 
unsteady cascade aerodynamic theory. 

Four external datasets must be furnished to ASTROP2 
for an analysis condition. These datasets define: (1) the 
steady aerodynamic inflow into each rotor, (2) the blade 
steady state geometry, (3) the blade dynamic 
characteristics in the form of mode shapes and frequencies, 
and (4) the blade setting angle and freestream 
aerodynamic information. In the present analysis, the 
aerodynamic solver of Ref. 9 was used to provide the 
nonuniform steady inflow (item(l)) into each rotor, and 
MSC/NASTRAN (Ref. 10) was used to provide items (2) 
and (3) as described below. Item (4) was user input to 
ASTROP2 directly. 

ASTROP2 combines the normal modes calculated by 
NASTRAN and a 2D-aerodynamic strip representation of 
the blade to calculate the generalized aerodynamic matrix, 
[A], described in Appendix A. It uses a frequency domain 
solution which predicts frequency, damping, and phase 
angle of aeroelastic modes. It is a single rotor analysis, 
meaning no rotor-to-rotor unsteady aerodynamics or 
structural dynamic interaction is accounted for. However, 
steady aerodynamic interaction between rotors can be 
accounted for by user input of the steady fiowfield into 
each rotor. This is the state of the art technology for 
aeroelastic analysis of CR propellers. 
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ASTR0P2 computes the unsteady aerodynamic loads 
(using the theory presented in Refs. 11 and 12.) and the 
damping and frequency for each respective interblade 
phase angle of the aeroelastically coupled modes. 

Code Modifications 

Previously, the ASTROP2 code modeled single 
rotation propfans with a uniform steady inflow velocity 
distribution, and the analysis was manually iterated to 
converge upon the least stable mode. Additions and 
modifications have been made to model each rotor of 
counterrotation propfans (neglecting structural and 
unsteady aerodynamic rotor-to-rotor interaction) and to 
automatically converge upon the least stable mode. The 
finite element geometry and corresponding modal 
deflection information for the aft blades are transformed to 
simulate the forward blades. In addition, the axial and 
circumferential nonuniform steady flow into each rotor 
can be input to include the rotor-to-rotor steady 
aerodynamic interaction effects. This nonuniform inflow 
velocity information is used to calculate the local relative 
Mach number for each aerodynamic strip and the 
unsteady aerodynamic loads. For a counterrotating 
configuration, the local Mach number distribution changes 
significantly due to aerodynamic interactions between the 
rotors and affects the unsteady aerodynamic loads. For 
single rotation configurations it is sufficient to assume 
uniform inflow velocity. A separate steady aerodynamic 
solver (Ref. 9) was used to provide the nonuniform steady 
inflow velocity information in the present analysis. For the 
blade designs analyzed in this paper, the use of a 
nonuniform inflow velocity distribution provided less 
conservative stability results than would have been 
obtained from a uniform inflow velocity. 

The aeroelastic stability analysis was modified to 
automatically iterate (to converge) and identify the least 
stable (damped) mode and the corresponding interblade 
phase angle for the rotor being analyzed. The procedure 
used is as follows. First, a value for blade vibration 
frequency is specified by the user. Using this value, 
ASTROP2 calculates the unsteady aerodynamic loads and 
aeroelastic eigenvalues for all the interblade phase angle 
modes. Based on the damping, the least stable mode is then 
identified. The frequency associated with this mode is then 
used to calculate new unsteady aerodynamic loads and 
aeroelastic eigenvalues. This procedure is repeated 
automatically till the frequency used to calculate the 
unsteady aerodynamic loads and the frequency associated 
with the least damped aeroelastic mode are the same. This 
routine saves a significant amount of time compared to the 
former manual method used. This is especially true when 
the relative tip Mach number exceeds 1 and convergence is 
slower. Also, the automatic iteration is useful in analyzing 
a large number of operating conditions as was required in 
the present study. 


Formerly, a Newton-Cotes method was used to 
integrate the unsteady aerodynamic forces on span wise 
strips to calculate the total lift and moment on the blade. 
The strip layout scheme in ASTROP2 lays strips normal to 
the blade leading edge. For propfan blades with varying 
sweep along the span, this scheme generates strips of 
unequal spanwise width. However, the Newton-Cotes 
method is only valid for strips of equal width and large 
errors are introduced when analyzing swept propfan 
blades with strips of unequal spanwise width. Therefore, 
an alternative integration method was needed. For 
increased accuracy in the ASTROP2 code, a trapezoidal 
integration scheme has been implemented to calculate the 
unsteady aerodynamic forces on blades. The trapezoidal 
integration method is valid for strips of unequal width and 
gives accurate results for swept propfan blades. Analyses 
using the modified code have been verified with existing 
propfan model wind tunnel flutter data from Ref. 1. 

Structural Model and Blade Geometry 

Blade geometry and material information are first 
input to a composite material and geometric preprocessor 
for NASTRAN, called COBSTRAN (Ref. 13). COBSTRAN 
produces grid mesh, element geometry, and material 
properties in the form of NASTRAN bulk data cards. 
MSC/NASTRAN then provides both: a geometric non- 
linear structural analysis to calculate the blade steady-state 
deformed configuration and total differential stiffness; and 
a complete vibration analysis which calculates the natural 
frequencies and mode shapes of the blade deformed state. 

Blade Description and Analysis Conditions 

The propfan models were designed for tests in the 
NASA Ames 14 ft. wind tunnel. The two blade sets under 
analysis were the primary blades for the program. The 
CMl blade set represents a geared propeller design and the 
CM2 blade set an ungeared propeller design. CMl blades 
are of higher aspect ratio than the CM2 blades and the 
design tip speed of the CMl blades is lower than that of the 
CM2 blades (see Table 1). The hub diameter of the rotors is 
0.216 m (8.5 in.). The design total thrust for the model 
blades at cruise RPM (two counterrotating blade rows, six 
blades per row) is 622.7 N (140 lbs). All the blades are 
made from layered pre-impregnated graphite composite 
material 0.08 mm (0.0032 in) thick, with 90 plies in each 
CMl blade and 82 plies in each CM2 blade. 

The finite element models for the CMl and CM2 
forward and aft blades are shown in Figures 2 and 3. A 
rigidly clamped connection was assumed at the blade 
shank end. The hub that retained each blade set was 
assumed rigid and was not modeled. 

The blade flutter analyses conditions for the CMl and 
CM2 blades are shown in Table 2. For both blade designs 
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there were three tip speeds analyzed at each of three Mach at 2.9%. For the aft blades, seen in Figure 5(b), the first 

numbers. In addition, the CM2 blades were analyzed at a mode shows the maximum percent stiffening at 15% and 

fourth Mach number at some conditions. Also, the analysis the third mode shows the minimum at 1.8%. 
at each operating condition was done for the first six 

natural modes of each of the four blade designs. Hence, Also shown in Figures 4 and 5 are the ranges of bench 

more than 216 conditions were analyzed for each of the measured natural frequencies of the manufactured blades, 

final blade designs. During the design process three other For both the forward and aft blades, there were 26 blades 

CM1 and four other CM2 designs were considered. Most measured using laser speckle interferometry. It can be 
of these designs were also analyzed over the full set of seen that the measured natural frequencies are typically 
analyses conditions. lower than calculated. The effect of these lower blade 

natural frequencies on calculated damping will be 
Results and Discussion discussed in the next section. 


The calculated results for the CMl and CM2 blade sets 
are presented in two parts: (1) the in vacuum blade natural 
frequency variation with rotational speed, and the mode 
shapes at the design rotational speed; (2) calculated 
aeroelastic damping plots , damping versus axial Mach 
number. In part (1) above, the calculations were made 
using MSC/NASTRAN. The calculations of blade steady 
state deflected positions included centrifugal loads but not 
steady airloads. For both the CMl and CM2 blade designs, 
the steady airloads were significantly smaller than the 
centrifugal loads. It was shown in Ref. 1 that the steady 
airloads do not have a significant effect on the natural 
frequencies. The real and imaginary parts of the 
eigenvalue for part (2) above were calculated using 
ASTROP2. Usually, 0.5% or greater structural damping is 
present in composite blades, but zero structural damping 
was assumed in the present analysis. This assumption 
makes the calculated stability results conservative. Six 
natural modes for each blade were included in the analysis 
although only results for the three least stable aeroelastic 
modes are presented here. 

Blade Natural Frequencies and Mode Shapes 

Figures 4 and 5 show the calculated variation of the 
first three natural frequencies with increasing rotational 
speed in vacuum for the CMl and CM2 forward and aft 
blades. The modal frequencies increase with an increase in 
rotational speed. This trend is expected because of 
centrifugal stiffening. 

For the CMl blades, the percent stiffening from 0 m/s 
blade tip speed to 228.6 m/s (750ft/s) is calculated for the 
modes shown in Figure 4. For the forward blades, seen in 
Figure 4(a), the first mode shows the maximum percent 
stiffening at 10.8% and the third mode shows the minimum 
at 2.5%. For the aft blades, seen in Figure 4(b), the first 
mode shows the maximum percent stiffening at 12.3% and 
the third mode shows the minimum at 3.9%. 

For the CM2 blades, the percent stiffening from 0 m/s 
blade tip speed to 274.3 m/s (900 ft/s) is calculated for the 
modes shown in Figure 5. For the forward blades, seen in 
Figure 5(a), the first mode shows the maximum percent 
stiffening at 15.5% and the third mode shows the minimum 


Figures 6 and 7 show the calculated iso-displacement 
contours of the CMl and CM2 blades for the first six 
natural modes at the design rotational speed. Note that all 
the modes of CMl and CM2 consist of coupled bending, 
torsion and edgewise motions. The modes referred to 
below are meant to identify the primary motion that can be 
associated with a single-degree-of-freedom beam mode. 
For both CMl blades, the first two modes are first and 
second bending respectively, and the third mode is torsion. 
Whereas, for both CM2 blades, the first mode is first 
bending, the second mode is first torsion, and the higher 
modes are platelike. 

Calculated Damning 

For a six bladed rotor, it is theoretically possible for the 
rotor to go unstable in six interblade phase angle modes for 
each of the six normal modes included in the analysis. 
Only the most critical interblade phase angle mode (the one 
closest to instability) is shown here for each of the three 
least stable normal modes. The unsteady aerodynamic 
loads were calculated using twenty chord wise strips across 
the blade span. 

Figure 8 shows the variation of the aeroelastic 
damping coefficient (the real part of the eigenvalue) with 
freestream Mach number for the CMl forward blades at 
three different blade tip speeds. The three least stable 
modes are shown. Figure 8(c) is discussed here in detail as 
it corresponds to the maximum tip speed analyzed for the 
CMl forward blade. In Figure 8(c), the interblade phase 
angles corresponding to the least stable modes are: first, 
180°; second, 60°; and fourth, 300°. A blade was considered 
to be flutter free when the aerodynamic damping was 
more than 0.002 in all the modes. From Figure 8, it can be 
seen that over the entire operating range, the CMl forward 
blades are flutter free. 

Figure 9 shows the variation of the aeroelastic 
damping coefficient with freestream Mach number for the 
CMl aft blades at three different blade tip speeds. The 
three least stable modes are shown. For the maximum tip 
speed analyzed, seen in Figure 9(c), the interblade phase 
angles corresponding to the least stable modes are: first, 
180°; second, 0°; and fourth, 300°. From Figure 9, it can be 
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seen that over the entire operating range, the CM1 aft 
blades are flutter free. 

Figure 10 shows the variation of the aeroelastic 
damping coefficient with freestream Mach number for the 
CM2 forward blades at three different blade tip speeds. 
The three least stable modes are shown. For the maximum 
tip speed analyzed, seen in Figure 10(C), the interblade 
phase angles corresponding to the least stable modes are: 
second, 120°; fifth, 300°; and sixth, 180°. Though the 
damping in the sixth mode is very small, it is considered 
aeroelastically stable due to the very high frequency. From 
Figure 10, it can be seen that over the entire operating 
range, the CM2 forward blades are flutter free. 

Figure 11 shows the variation of the aeroelastic 
damping with freestream Mach number for the CM2 aft 
blades at three different blade tip speeds. The three least 
stable modes are shown. For the maximum tip speed 
analyzed, seen in Figure 11(c), the interblade phase angles 
corresponding to the least stable modes are: fourth, 60°; 
fifth, 300®; and sixth, 240®. From Figure 11, it can be seen 
that over the entire operating range, the CM2 aft blades are 
flutter free. 

The results presented in Figures 8-11 were obtained 
using the calculated blade natural frequencies. As seen in 
Figures 4 and 5, the measured blade frequencies were 
typically lower than the calculated frequency values. To 
study the effect of these lower frequencies, the analyses 
were repeated for natural frequencies 10% lower than the 
calculated ones. The aerodynamic damping remained the 
same for most of the modes and increased in some cases. 
All the blade designs were stilt predicted to be flutter free 
within the planned operating range using these lower 
natural frequencies. However, although aerodynamic 
damping was not changed much in these blade designs by 
the lower natural frequencies, in other designs 
aerodynamic damping may be adversely affected. Hence, 
each design must be independently checked. 

The aeroelastic modal frequencies were very close to 
the natural frequencies for all these blade designs. This 
indicates that there was very little aeroelastic coupling 
between different modes and all the aeroelastic modes 
were well separated. 

In the wind tunnel test, the CM1 and CM2 blade sets 
did not flutter. This enhances the confidence in the 
ASTROP2 code as an aeroelastic design tool for both single 
and counter rotation propfan configurations. 

Conclusions 

An aeroelastic stability research code for single 
rotation propfans was modified and successfully applied to 
counterrotation propfans. 


The code was used to analyze the aeroelastic stability 
of two different counterrotation propfan blade designs for 
the joint Navy/Air Force/NASA cruise missile wind tunnel 
test. Both the blade designs, called CM1 and CM2, were 
predicted to be flutter free. Experimentally this was 
confirmed over the entire operating range of the wind 
tunnel test . 
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The generalized aerodynamic force vector can be 
expressed in terms of lift, moment, normal modes, and 
normal coordinates as 


{Q) = 7tpco : 


H 


L 

M a 


ds =[A] {q> 


(A4) 


Appendix A 

Calculation of Generalized Aerodynamic Forces 


It is assumed that the lift and moment per unit span 
due to blade motion are linearly related to the 
displacements and their derivatives with respect to arc 
length along a reference line. The blade is divided into a 
series of discrete aerodynamic strips. Each strip has two 
motions, pitching (a) and plunging (h) about an arbitrary 
reference line. The lift and moment expressions per unit 
span are 



= Tip CO 2 


111 1 12 - ll6 
-l2ll22 - ^26. 


h 

a 

h’ 

a' 

h" 


a" 


(Al) 


The prime denotes the differentiation with respect to the 
arc length along the reference line. The elements l n , 1 12 , 
etc, are defined in Ref. 1. The pitching and plunging 
displacements and their derivatives are given by 


where 


h 
a 
h' 
a' 
h" 
\ a” 


I [w] 

= [W] 
([W-] 


<q> 


[W ]= w, w 2 w m 

- Al A 2 A™ 


(A2) 


(A3) 


where 1 represents the length of the reference axis. 
Substituting Eqs. (Al) and (A2) into Eq. (A4), 


[A] = Tip (JO 2 


[W] T [LCOEF] 


Jo 


/ [W] \ 
[W] 
[wl 


ds 


where 


[LCOEF] = 


I 11 1 12 Le 
_ I2 1 122 — ^26. 


(A5) 


(A6) 



CMl 

CM2 

Freestream Mach 

0.4 to 0.85 

0.4 to 0.9 

number 



Max. Rotor Speeds 

10417 

14475 

rpm 



Tip Diameter 

0.425(16.75) 

0.368 (14.5) 

m (in) 



Design Conditions 



Mach number 

0.7 

0.7 

tip speed, m/s (ft/s) 

213.4(700) 

274.3(900) 

rpm 

9580 

14225 

blade angle, 3/4R, 0 



. forward blade 

58.5 

51.8 

aft blade 

56.4 

50.0 


TABLE 1 - BLADE DESIGN CONDITIONS 


Blade 

Mach Numbers 

Tip Speeds, m/s (ft/s) 

CMl Forward and CMl Aft 

0.4, 0.7,0.85 

137.2 (450), 182.9 (600), 228.6 (750) 

CM2 Forward and CM2 Aft 

0.4, 0.6, 0.8, 0.9 

182.9 (600), 228.6 (750), 274.3 (900) 


TABLE 2 - BLADE FLUTTER ANALYSES CONDITIONS 
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10417 RPM 





NATURAL FREQUENCY (Hz) NATURAL FREQUENCY (Hz) 











Aerodynamic Damping Aerodynamic Damping Aerodynamic Damping 



(a) 



(b) 



(c) 

Figure 8 Variation of damping with Mach number and blade 
tip speed for CM1 forward blade (a) 137.2 m/s (450 ft/s) 
(b) 182.9 m/s (600 ft/s) (c) 228.6 m/s (750 ft/s) 



(a) 



(b) 



(c) 

Figure 9 Variation of damping with Mach number and blade 
tip speed for CM1 aft blade (a) 137.2 m/s (450 ft/s) 
(b) 182.9 m/s (600 ft/s) (c) 228.6 m/s (750 ft/s) 












Aerodynamic Damping Aerodynamic Damping Aerodynamic Damping 



(a) 



(a) 




(b) 


(b) 



(c) 

Figure 10 Variation of damping with Mach number and blade 
tip speed for CM2 forward blade (a) 182.9 m/s (600 ft/s) 
(b) 228.6 m/s (750 ft/s) (c) 274.3 m/s (900 ft/s) 



(c) 


Figure 11 Variation of damping with Mach number and blade 
tip speed for CM2 aft blade (a) 182.9 m/s (600 ft/s) 
(b) 228.6 m/s (750 ft/s) (c) 274.3 m/s (900 ft/s) 
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